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reduced, but the break-up of the jet occurs further from the orifice, and for 
the lowest pressure not until after the lip of the oscillator or pipe has been 
passed. In this case, as was before mentioned, only a small amplitude of 
vibration is maintained, and the effect of the dynamic stiffness of the jet 
preponderates. 

Thus the fact that low wind or water pressure depresses the pitch, while 
higher pressures raise it, is accounted for, although it is impossible to calcu- 
late the amount from purely dynamical data (viz., velocity, pressure, density, 
and viscosity), as no theory is at present capable of dealing with the formation 
of eddies or the final result of fluid instability. 
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The calculation of the whirling speeds of loaded shafts of turbines and 
turbine-driven machinery is usually rendered complex by reason of the 
varying distribution of the load along the length of the shaft, as well as by 
variations from place to place in the shaft's diameter. 

Of the several methods in use for the estimation of the first whirling speed 
in such cases, probably the best known is that based on the Dunkerley 
hypothesis. This method, which was put forward by Dunkerley in 1893* 
and illustrated by him in many applications, consists in dividing the load 
into a convenient number of parts acting at appropriate positions along the 
shaft, and determining by means of simple formulae the whirling speed for 
each such part of the mass regarded as acting alone. It is then assumed 
that the sum of the squares of the reciprocals of these partial whirling speeds 
is equal to the square of the reciprocal of the first whirling speed for the 
whole load. 

Expressed algebraically, if g>i, 0)2, ws, etc., be the whirling speeds of the 
shaft loaded in turn with the separate masses regarded as acting alone, and 

^ <PhiI. Trans.,' A, vol. 185, pp. 279-360 (1894), especially §§28, 61, 62. 
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XI be the first whirling speed of the loaded shaft, then Dunkerley's empirical 
rule consists in the assumption that 

Q? G>1^ G>2^ (Ji)^ 

In many cases this method gives excellent results, and, having regard to 
the uncertainty of the data and the purposes for which the calculation is 
required, quite sufficient accuracy is often obtained. But the method is not 
invariably safe, and in certain applications may give results very seriously in 
error. 

That the empirical rule is subject to such limitations was known to 
Dunkerley, who, in the paper referred to, proposed, in certain cases in which 
the calculated result did not agree with experiment, to substitute a formula 
of the type (for two loads) 

1 1 . a 



12^ ft)i^ ft)2^' 

where a has a value less than unity.* 

The need for such modification in certain cases was also pointed out by 
Chree.f It is not, however, clear beforehand what cases require modification, 
or the extent of the modification, i,e, the value to give a. 

It is proposed in the following note to show how Dunkerley's empirical 
rule may be derived rationally, and as a consequence to show its limitations, 
and how it may be extended to give more accurate results in any application. 

It is to be noted that the method used below is of general application, but 
to fix the ideas and render the argument clearer, the case will be taken of a 
uniform and light shaft freely supported by bearings at its ends, and loaded 
with a number of discrete masses mi, m2, ..., m„, of negligible moment of 
inertia at definite points along its length. 

As above, let coi, 0)2, ..., ^w, be the whirling speeds corresponding to each of 
the masses mi, m2, ..., m„, acting alone ; and let Oi, XI2, ..., O^, be the first, 
second, etc., whirling speeds of the loaded shaft with all loads in position, 
corresponding to the different possible types of vibration. In any assumed 
possible type of lateral bending of the shaft under the action of centrifugal 
forces and subject to the end restraints, let iii, U2, ..., '^^m, be the lateral 
displacements from the straight form at the points where the masses mi, 
^2, ..., Wm, are fixed. 

In the first instance, consider the shaft subjected simply to a force Xi 
acting at right angles to its length Z at a point distant Zi from one end. Let 

* Dunkerley, ibid., §§ 30, 33, 44. 
t ' Phil. Mag.,' May, 1904, p. 641. 
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E be Young's modulus for the shaft and I its moment of inertia of cross- 
section about a diameter. Then, following the BernouUi-Eulerian method, 
since u is the displacement at the point z parallel to X, we have, for the 
shaft from ;<; == to 2; = %. 

«~EI^ = X'^, 

where. X^ is the force exerted at the bearing at = 0. 
Similarly for the portion of the shaft zi to I we have 

^■EltJL^'X.''{l-z% 

dz" 

where X" is the force exerted at the bearing at z = L 
Integrating we obtain 

QMu = 6EI% ^ + X'^ (zi^-z^) 

Zi 

and • 6ELu' = 6llt^i j^^ + Xf' (l-z) (z^zi) (21^ 

i — Z\ 

where % is the value of 1^ at zi, Now, neglecting "th- 
shearing," the curves have a common tangent at Zi, Henc 

But X = Xi — J — , and X = Xi y. 

Hence 6EIte = Xi2j(l—%)(2fei— %2— #), 

and %WLlw' ^Xizi(l-z){2lz-z^-zi^), 

These equations give the displacement at any point along the shaft conse- 
quent on the action of the force Xi at 01. 

If now a number {n) of such forces X act at several points along the shaft, 
the deflection at any point will be the sum of the separate defiections at that 
]point due to each of the forces acting separately. Accordingly we obtain 



8 = 111 



6EK%= S XsZs{l—Zr){2kr'-Zr^'-'Z^^)-^ 2 X^Zril — Zs) {^Iz^—Zs^ --Z,^), 

• where the first summation includes all the forces on the shaft from to Zr, 
and the second summation includes all those from ^^ to ^. 

If, now, the forces X^, etc., be due to the centrifugal action of the masses 
m^, etc., when the shaft is rotating with uniform angular velocity XI in a 
deflected form, %, etc., being the corresponding displacements, then 

Xs = m^sO^i^s, etc. 

If the shaft is not rotating, when determining the speed of lateral 
vibration we put X^ = mjD?iis in the same way, where X2 is now the speed of 
lateral vibration. 
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Substituting we obtain 



o 



s =.r 



s = 



s ~ fb 



+ 0^ % msZr(l—Zs)(2ks—Zs^--Zr^)Usf 

where the first summation extends to each of the masses up to and including 
nir towards one end of the shaft, and the second summation to those masses 
beyond 7% towards the other end. 



Putting 



6EK 



#, the above equation may be written in the more 



convenient form 

where the coefficients /i,/2, etc., are functions of the magnitudes and positions 
of the masses mi, mg, etc., and of Zr, 

Taking m^ to be each of the masses mi, m2, ..., m„, in turn, we have the 
equations 

^Ui = aitCi + a2tl2 + . . . + Ctntln, 
^U2 = hUi + h2U2 + . . . + bnUn, 



^qi^ = hiUi + ^2^2 + . . . + hnUn. 



It is worthy of note that the coefficients are related to one another in a 
symmetrical manner, viz. 

^ __ m2 cts __ ^3 



-i = ~i, etc., 

C2 m2 



Ml Ci mi 

a result to be expected from the well known theorem that the displacement 
at a point A of a beam due to a load at another point B, is equal to the 
displacement at B due to the same load acting at A.* 
Eliminating Ui, U2, ..., Un, we obtain the determinant 

Cti — ^, (X>2) .♦. Cl-n 

h, h—^, ... K 



= 0. 



The roots of this equation in # give the n whirling speeds of the shaft 
loaded with the n masses mi, m2, ..., m„, all on together, viz., Oi, O2, ..., f2„. 

This is not the ordinary determinant which gives the periods of vibration 
of a connected system of masses. The present determinant involves 1/W 
instead of O^ as in the ordinary one, from which the same result may be 

* Eayleigh, * Collected Papers,' vol. 5, p. 170. 
VOL. XCV. — A.. K 
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obtained. The present form, however, is more convenient for the numerical 
computation of whirling speeds. 

Expanding the determinant in powers of $ we obtain 

Thus we find for the sum of the roots 

If k = 6EI^, we have 

Next let us consider each mass acting alone. We find, either directly, or 

from the above general equations on equating each of the remaining masses 

to zero, 

(fyih = cciUi, (f>%i2 — h2U2, etc. 

Hence (jbi = ai, <f)2 = &2, etc. 

Also (f)i = — 5, 62 = ~„ etc. 

We thus see that 



(Dn 






03 - ^ ^• 



Thus the sum of the squares of the reciprocals of the several whirling 
speeds of a shaft carrying a number of loads i§ equal to the sum of the 
squares of the reciprocals of the whirling speed of each mass regarded as 
acting alone. 

Now, Oi, ^2^ f^3> ...J ^n, being the several whirling speeds of the loaded 
shaft, are an increasing series of numerical values. If O2, O3, etc., are 
considerably greater than Oi, we may to a first approximation neglect the 
sum of 1/122^ 1/^3^ etc., compared with l/Oi^. We then have, approxi- 
mately, provided we are justified in the omission of O2, O3, etc., 

1__ 1 , _1_ , , 1 



Oi^ w^ mi *" o)/ 

This is Dunkerley's empirical rule. 

Its degree of approximation then clearly depends on the values of the 
second and higher whirling speeds of the loaded shaft relatively to the first. 
Necessarily, it gives too low a value for the first whirling speed. 

To obtain a numerical idea of the degree of approximation, it will be 
convenient to consider the case of a uniformly loaded and uniform shaft, in 
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which the equations giving the several whirling speeds under different 
conditions of support are well known. 

Let g stand for the ratio -tt-ttt-s-^- 

(l/Oi^) 

A. Taking first the case of a shaft fixed in direction at one end and free at 
the other. Then* the natural frequencies must satisfy 

cos />t cosh yU. "+ 1 = 0, 

where in our application fxP is proportional to 12. 

Accordingly, the frequencies of the tones or the values of the several 
whirling speeds are proportional to 1, 6*27, 17'55, 34*39, etc. 

rurther,f the expansion of the equation in fx leads to 

1— jfi.4. = 

1 1 

SO that S — r = 



/^* 12 
But //,* is (1-875104)4 or 12-3624. 

Hence q = ^^^ = 1*0302. 

B. Consider next the case in which the shaft is simply supported at each 
end. Then| the natural frequencies must satisfy sin yu, = 0. 

Thus the frequencies of the tones or the values of the several whirling 
speeds are proportional to 1, 4, 9, 16, etc. 

Also in this case the expansion of the equation in fx is 

so that t --■^, t —2—2 = "^ ; therefore t — = jr^. 

fjb^ 6 //'I yU'2 120 )Lt* 90 

But jJbi = TT. 

Therefore q ^ ^^ V ^? = ^ = 1*0^23. 

^ (l//.i4> 90 

0. If the shaft be fixed in direction at one end and supported at the other, 

the equation for jx is§ 

cot jjb = coth fi, 

which on expansion gives 1— ^^+... = 0. The frequencies are pro- 
portional to 1, 3*24, 6*76, 11*56, etc. 

■^ See Rayleigh, 'Theory of Sound,' vol. 1, § 173, and Dunkerley, ihid.^ §9. 

t Eayleigh, ihid.^ § 175. 

t Eayleigh, ihid.^ §181 ; Dunkerley, ihid,^ §10. 

§ Eayleigh, ihid.^ §§ 180, 174 ; Dunkerley, ihid.^ §14. 
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Also S —7= ^777: and jii = 3*9266, 

fjL^ 210 

so that /i/ = 237-721. 

Therefore a = ^^^- = 1*132. 

D. Shaft fixed in direction at each end. 

The equation for yu, is* 

cos /x, cosh /^—l = 0, 

a* 
which reduces to 1 •— 7^ + . . . = 0. 

420 

The frequencies are proportional to 1, 2*76, 5*40, 8*93, etc. 
ence 2^ ■— : = 



4 420 



But 1X1^ = (4-73004)*t = 500-564. 

Therefore a = ^,9'J ^^} = I'l^^. 

E. Finally take the case of a shaft simply supported in three bearings, the 
spans being of equal length. 

Then| fju must satisfy either 

sin /i = 0, or cot //< = coth fi. 

The frequencies are proportional to 1, 1*56, 4, 5*06, etc. In the former 
case 2(1/m*)= 1/90, as in (B) above; for the latter :S(l//i*) = 1/210, as 
in (G). For both sets of frequencies, therefore, we have %(1//jl^) = 1/63. 

But fJii = TT. 

Hence q - -^^^^ ^ 63 ^^ 

Summing up these results, we see that the values of Oi obtained by 

omitting the higher whirling speeds are thus too low, in the following 
ratios * ' '■ ■ 

0-985, 0-961, 0-940, 0-916, 0'804. 

Thus in some of these cases the agreement is good ; in the last it is quite 
insufficient, and this is not necessarily the worst case that may arise. 

These results may also be obtained otherwise, since we know for these 
several systems of support the whirling speed of a massless shaft carrying an 
isolated mass at any point along its length. A uniformly loaded shaft may be 

"^ Eayleigh, ibid., §172 ; Dunkerley, ibid.^ §17. 
t Kayleigh, ibid., §174. 
I Dunkerley, ibid., §15. 
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considered made up of an infinite number of elements of mass spread along 
the shaft, and we can, therefore, by integration obtain the value of ^(l/co^), and 
compare it with the exact value of l/Oi^ in each case. 

To obtain a better approximation in cases where O2 is not much greater 
than ^Ql, and where in consequence the Dunkerley rule is not sufficiently 
exact, the preceding determinant is to be expanded to a further power of <l> 
thus 

$w-.<j)«-i(ai + 52 + C3+...-f A^) 

+ ^"""^ {{ciih2) + (aiCs) + . . . + (hcs) + . . . + {g7i~ 1 Z^??)} — • . . = 0, 
where (ai&2) stands for the determinant 

As before, the sum of the roots #1 + $2 + . . . + ^j? is equal to ai 4- 62 4- ... 4- ^hi- 
Also the sum of the products of the roots taken two at a time 

^1<J>2 4- <I>1^3 4- . . . 4- $2<E>3 + . . . + #w- l*n 

= (aih) 4- (aiCs) 4- ... 4- (^2^3) 4- ... 4- {On-ihi)- 

7 7 

But ^n = TTgj ^^^d hn = — 2 as before. 



■"n ^n 



Also (^162) is the product of the first and second whirling speeds of mi and m2 
acting together, but with the other masses omitted ; and so on. 
We thus find that 

S — K = S -1^ and 2 ^ ^^ ^ = S 



where «', &>" are the first and second whirling speeds when the masses are 
taken as acting two at a time. 

If now we may neglect XI3 and the higher whirling speeds we have 

1 4-_Ji = 2 -^ 



X ^ -L 

— . ^ 



and from these two equations a more exact value of Oi may be deduced, 
together with an approximate value of O2. 

In this case it becomes necessary to determine the whirling speeds for each 
of the masses regarded as acting alone, and also the first and second whirling 
speeds for the masses acting two at a time in all possible pairs. 

Dunkerley's empirical method may thus be extended to the following 
rule : — 

Find the whirling speed (00) for each load acting separately ; and the 
first and second whirling speeds {(o\ &)'') for each combination of two loads 
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acting without any of the others. Then if Oi and O2 are the first and second 
whirling speeds for the whole system, they are given by 

1 I 1 I_- V 1 



Oi^ * Os^ • (1)2' 



If we expand the determinantal equation to further powers of ^ we will 
find in the same manner that the sum of the squares of the reciprocals of the 
products of the n whirling speeds taken r at a time is equal to the sum of the 
squares of the reciprocals of the products of the r speeds for each set of r 
masses, the other n—r masses being supposed reduced to zero. 

Further, as a similar determinantal equation will obviously be obtained in 
any general case in which linear displacement equations of the preceding 
type apply, we may express with greater generality the results arrived at 

Thus the sum of the squares of the reciprocals of the products of the n 
speeds of vibration taken r at a time of n masses connected elastically, is 
equal to the sum of the squares of the reciprocals of the products of the 
r speeds of vibration for each set of r masses, the other ^— r masses being 
supposed to be obliterated. 

The particular case when r = 1 is perhaps the most interesting. The sum 
of the squares of the reciprocals of the n speeds of vibration of the system of 
n masses is equal to the sum of the squares of the reciprocals of the speeds of 
vibration of each of the n masses regarded as acting alone. 

Further, if the n masses be divided into any number of groups the sum of 
the squares of the reciprocals of all the speeds of vibration of each of the 
groups, taken for all the groups, is equal to the sum of the squares of the 
reciprocals of the n speeds of vibration for all the n masses acting together. 
The theorem admits of other generalisations. 

In general, then, the n whirling speeds are the roots of the determinantal 
equation obtained above. In any particular application numerical values 
may be substituted for the masses, their positions, and the elastic constants 
of the shaft. The roots of the resulting numerical determinantal equation 
may be evaluated by successive approximation. Usually not more than the 
first two or possibly three speeds are required. 

This determinantal method, although tedious when the number of isolated 
masses is great, is straightforward and exact. It may be made to take 
account of variations of shaft diameter, and also of the moments of inertia of 
the loads, and of different bearing conditions. It is quite general and makes 
no assumption of an approxin^ate displacement curve of gravitational or other 
type. 
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In cases where with sufficient accuracy whirling speeds higher than the 
second may be neglected we have, in the foregoing notation, 

11. 

4- -TT-^ ==7 C)^i + &2 + ^3 + • • • + hn, 



^^^ TTTFT^ ^ (^^1^2) + {cilC^) + . . . + (&2^3) + . . . + {gn- 1 K) ] 



am 



2 



i.e, the sum of all the second order determinants about the principal 
diagonal. 

Oi and XI2 ri^ay thence be determined. 
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During the last few years some attention has been paid to the spectrum of 
gases in the neighbourhood of the cathode discharge, and, therefore, it is 
thought to be worth while putting on record some experiments that were 
made about ten years ago, but which were not considered at the time to be of 
sufficient importance to be published, as they led to no definite conclusion. 
One of the authors* had noticed that, in a vacuum tube containing helium 
and mercury vapour, the orange mercury line 6152 appeared very brightly, 
although it is not usually seen at all in a tube containing a diatomic gas. 
Continuing this line of research, other metals were experimented on, with 
those metals as electrodes. It was at once noted that in many cases, at 
comparatively high pressures and *low current densities, lines in the 
spectrum of the metals composing the electrodes were visible in the 
neighbourhood of the cathode, and often in other parts of the tube as well. 

The following experiments were made, with the view of studying more 
thoroughly this effect. Cadmium was the metal chosen, on account of the 
ease with which the metallic spectrum was produced, and also because the 
cadmium bright lines extend over a considerable portion of the visible 
spectrum. Other metals that behave in a similar manner are zinc, 
magnesium, lead, bismuth. The gases experimented with were helium, 
neon, argon, and xenon. 

^ ' Boy. Soc. iProc.,' vol. 71, p. 25. 



